IPPP/02/83 
DCPT/02/166 



From weak-scale observables to leptogenesis 



X 



Sacha Davidson 

IPPP, University of Durham, Durham, DHl 3LE, UK 



Abstract 



m 
o 
o 

■ Thermal leptogenesis is an attractive mechanism for generating the baryon asymmetry of 

the Universe. However, in supersymmetric models, the parameter space is severely restricted 
by the gravitino bound on the reheat temperature Tj^h. For hierarchical light neutrino 
' masses, it is shown that thermal leptogenesis can work when T^h ~ 10^ GeV. The low- 

\ energy observable consequences of this scenario are BR{t €7) ^ 10^^ — 10^^. For higher 

, TnH, thermal leptogenesis works in a larger area of parameter space, whose observable 

' consequences are more ambiguous. A parametrisation of the seesaw in terms of weak-scale 

_i I inputs is used, so the results are independent of the texture chosen for the GUT-scale Yukawa 

^ ■ matrices. 
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2 ! 1 Introduction 

i-S^ ' Leptogenesis is an appealing mechanism for producing the baryon asymmetry of the Universe 151. In 

the seesaw model|3], heavy singlet ("right-handed") neutrinos vn decay out-of-equilibrium, producing a 
^ 1' net lepton asymmetry, which is reprocessed by Standard Model (SM) B + L violating processes 0] into 

a baryon asymmetry. A natural and cosmology-independent way to produce the vn is by scattering in 
the thermal plasma. This scenario is refered to as "thermal leptogenesis" . However, the lightest can 
be produced only if their mass Mi is less than the reheat temperature T^h of the plasma after inflation. 
In addition, the asymmetry is proportional to Mi (5], so there is a lower bound on Mi to get a large 
enough asymmetry. This implies 10^ GeV < Mi < Trh- 
d ' The seesaw is an attractive minimal extension of the SM that generates the observed small v masses. 

Three right-handed neutrinos, with large majorana masses Mi, are added to the Standard Model, along 
with a Yukawa matrix for the neutrinos. It is desirable to supersymmetrise the seesaw, to address the 
hierarchy between the weak scale and the Mi. In the SUSY seesaw, T^h must be low enough to avoid 
over-producing gravitinos |^[7] — the canonical bound for gravity mediated SUSY breaking is Tjiu ^ 10^ 
GeV. The aim of this paper is to identify the parameter space where thermal leptogenesis can work, 
taking Ml - Trh 10^ GeV. 

The SUSY seesaw has more low-energy consequences than the non-SUSY version. It induces lepton 
flavour violating (LFV) entries in the slepton mass matrix, which can lead to radiative lepton decays[S], 
such as fi 67, at experimentally accessible rates. Eighteen parameters are required to define the 
neutrino and sneutrino mass matrices (in the charged lepton mass eigenstate basis), which is the same 
number as there are high scale inputs for the seesaw model 9 . It can be shown that the SUSY seesaw can 
be parametrised with the sneutrino and light neutrino mass matrices in a texture model independent 
way. That is, the high-scale physical inputs of the SUSY seesaw — the masses Mi and Yukawa coupling 
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Yi, — can be "reconstructed" from the neutrino and sneutrino mass matrices^. The baryon asymmetry 
can therefore be expressed as a function of weak scale observables. In this paper we identify the 
ranges of experimentally measurable quantities which are consistent with thermal leptogenesis. This 
phenomenological analysis differs from previous work [TTllT^OiniTHliraiTTllT^iroil^iniO^ bv 

making minimal assumptions about the high scale theory: we assume the SUSY seesaw and universal 
soft masses at the GUT scale. The usual approach is to assume a GUT-scale texture that generates the 
desired neutrino mass matrix, and discuss leptogenesis — the aim here is to input the slepton mass matrix 
instead of a texture. This analysis should be consistent with all GUTS and texture choices covered by 
these assumptions. 

Section 2 includes notation, and a review of leptogenesis and our parametrisation of the seesaw 
model. Section 3 presents approximate analytic expressions for the quantities on which leptogenesis 
depends. The low energy signatures of the parameter space where thermal leptogenesis works are 
discussed in section 4. CP violation is briefly discussed in section 5. The results are discussed and 
summarised in sectional 

2 Review 

The observed deficits in muon neutrinos from the atmosphere | 23| and in electron neutrinos from the 
sun^4, 25 , 2^ can be fit with small neutrino mass differences. The recent KamLAND observation of a Ve 
deficit from reactors confirms the neutrino mass explanation of the solar neutrino puzzle P?]. The small 
Am^ are consistent with three patterns of neutrino mass: hierarchical (Am^^^ = m|, Arrig^; = vn^, 
degenerate (ma ~ m2 ~ m\ 3> Am^^^) and quasi-Dirac (m| ~ — Arn^^^, Am^^; = mg — m|). The 
leptogenesis scenario considered in this paper, where the vb. are produced by scattering in the plasma, 
does not work for degenerate |3] (see also |2H1 for a detailed discussion). The quasi-Dirac spectrum 
could be interesting, although it is possibly disfavoured by supernova data [^Hl- We assume the are 
hierarchical, so the neutrino masses are much smaller than the charged lepton and quark masses. These 
small masses can be naturally understood in the seesaw model. 

In subsection l2.1l the seesaw is reviewed from the top-down; introducing new physics at a high scale 
Mxi and seeing its low energy implications. This approach has been followed by many model builders 
who construct a natural or symmetry-motivated structure of the high-scale mass and Yukawa matrices, 
and then study its low energy consequences. See e.g. [21] for early works that produce neutrino mass 
matrices with small mixing angles, and |31) for more complete up-to-date references. Lepton flavour 
violation due to the SUSY seesaw, which could be observed in 1^ li^ 13122] or in slepton production 
and decay at colliders has also been extensively studied from a top down approach (see e.g citations 
of Recent studies (for instance [33]) have considered the branching ratios for £j — > Ici in models 

that induce the two observed large mixing angles among the light leptons.^. 

Subsection 12.31 is a "top-down" review of leptogenesis, which is the obvious approach [111 I12L [T51 
El El El El 1201 1^ EH since the asymmetry is generated at high scales. See |TT] for examples and 
models. The translation between this approach and our bottom-up phenomenological analysis is not 
obvious, so it is difficult to relate our work to these papers. A phenomenological analysis of leptogenesis 
in non-SUSY (so no LFV) SO(IO) models was discussed in JHIES, with particular attention to the 
low-energy CP violation. A Yukawa-matrix independent analysis has also been done in the case where 
there are only two right-handed neutrinos [T7|. 

^This "reconstruction" would require universal soft masses at the GUT scale, and improbable experimental accuracy 
at the weak scale, so is in practice impossible. 

■^see also |35| for a more phenomenological discussion of — ► 67 
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2.1 Notation and Numbers 



We consider the supersymmetric see-saw for two reasons: first, supersymmetry stabilizes the Higgs mass 
against the quadratic divergences that appear due to heavy particles {e.g. the right-handed neutrinos). 
Secondly, the slepton masses enter our bottom-up parametrisation of the see-saw. 
The leptonic part of the superpotential reads 

Wi^j, = e^^YeL • Hd + v'rY.L ■ - ^i^f^Miy^ (1) 

where Li and Cfn {i = e,fj,,T) are the left-handed lepton doublet and the right-handed charged-lepton 
singlet, respectively, and Hd {Hu) is the hypercharge —1/2 (-f 1/2) Higgs doublet. Yg and Y^, are the 
Yukawa couplings that give masses to the charged leptons and generate the neutrino Dirac mass, and 
is a 3 X 3 Majorana mass matrix. This is the minimal seesaw; additional terms are possible, for 
instance in SO(IO) models a small triplet vev (T) is probable j36j. leading to a vl{T)i^l mass term. 

We work in the left-handed basis where the charged lepton mass matrix is diagonal, and in a basis 
of right-handed neutrinos where is diagonal 

DM^diag{Mi,M2,M3), (2) 

with Mi > 0, and Mi < M2 < M^. In this basis, the neutrino Yukawa matrix must be non-diagonal, 
but can always be diagonalized by two unitary transformations: 

Y, = V^DyVl, (3) 

where Z?y„ = dmg(yi, j/2, J/s) and 2/1^2/2^ 2/3- Later in the paper, we will assume that Dy^ is 
hierarchical, with a steeper hierarchy than is in the light neutrino mass matrix: (2/1/2/2)^ ^ mi/m2. 

It is natural to assume that the overall scale of A4 is much larger than the electroweak scale or any 
soft mass. Therefore, at low energies the right-handed neutrinos are decoupled and the corresponding 
effective Lagrangian contains a Majorana mass term for the left-handed neutrinos: SCiep = —^v^rrii^vA- 
h.c, with 

= niD^X^'mD = YJM-^Y,{HIY. (4) 

We define the Higgs vev (-ff^)^ — ''^'u — sin^ /?, where u = 174 GeV. In the basis where the charged- 
lepton Yukawa matrix. Ye and the gauge interactions are diagonal, the [mj,] matrix can be diagonalized 
by the MNS i3S| matrix U according to 

U'^[m^]U ^ diag{mi, 7712,1713) = D„i^, (5) 

where /7 is a unitary matrix that relates flavour to mass eigenstates 

(6) 

and the rm can be chosen real and positive, and ordered such that toi < m2 < m^. Assuming hierarchical 
left-handed v masses, we take m§ = Am^^^ = 2.7 x lO^^eV^ ^ and ml = '^m^oiai = 7.0 x 10~^eV^ 
gni- This corresponds to = 5.2 x 10"^ eV (3.9 - 6.3 x lO'^ eV at 90% C.L.), and 7712 = 8.2 x lO'^ 
eV (7 — 15 X 10^"^ eV at 3 ct). mi is unknown, usually unimportant, and we take it to be 7712/10. As we 
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shall see, the baryon asymmetry is weakly dependent on tan/3 in the parametrisation we use, so we set 
sin/3 = 1. 

U can be written as 

C/ = y •diag(e-^"^/2,e-''^'/2,l) , (7) 
where (j> and (/>' are CP violating phases, and V has the form of the CKM matrix 

C13C12 C13S12 5136"**' 

V = { -C23S12 - 5235136126*'' 623C12 - 5235135126*'' 523C13 | • (8) 
523512 - 6235136126*'' -S23612 - 6235135126*'' 

The numerical values of the angles are .28 < tan^ 9soi < -91 (Scr), with best fit point tan'^ dgoi = -44 
|40|. so 6soi = -41. We take Oatm — tt/4:. The CHOOZ angle 6*13 is experimentally constrained sin0i3 < .2 
|41| . Considerable effort and thought has gone into designing experiments sensitive to smaller values 
of 013. J- PARC hopes to reach 0(0.05) ^42,, and a neutrino factory could detect ^13 as small as 

0.02-^ 0.001 1321 Eg. 

We assume a simple gravity-mediated SUSY breaking scenario, with universal soft masses at the 
scale Mx- The sneutrino mass matrix (in the charged lepton mass eigenstate basis) can be written in 
the leading log approximation as 

(ml)^^ ^ (diagonal part) ^+ " [Yt],,[Y.],, In ^ , (9) 

where "diagonal-part" includes the tree level soft mass matrix, the radiative corrections from gauge and 
charged lepton Yukawa interactions, and the mass contributions from F- and D-terms. 
The branching ratio for ij — > £i"f can be estimated 

BRje.^ia) - ~ ,2,...r,-7~2y. y , / lOOGe^ ^ / tan /3 " 

(10) 

where C ~ 0(0.001 0.01), and Vl diagonalises the second term of eqn 0. More accurate formulae for 
the branching ratios can be found in [32] . To further simplify these estimates, it would be convenient 
to assume that Vl — Vl- That is, the lepton asymmetry will be a function of the angles of Vl, and 
it would be simplest to estimate £j £^7 using the angles of Vl for those of Vl- This a reasonable 
approximation when Blij ^ ^^Rij < i)i where 9B,ij ( ^iij) is an angle of VniVL)- For hierarchical 
Yukawa eigenvalues, this is likely to be true, even if an angle in Vr is large, because the usual texture 
estimate for Olij would be ^Jyi/ yj . We assume this condition is verified, so the principle contribution 
to [ml],j is (xylVli^^VLkj- 

Table 12.11 lists the current upper limits on the /i ej, t — s- 67 and r branching ratios, 

and the corresponding bounds on Vlsj that can be estimated from eqn H10|l |44| . It also contains the 
hoped for sensitivity of some anticipated rare decay searches. Colliders could also be sensitive to flavour 
violating slepton masses pS] . 

Leptogenesis will depend on angles of Vl- In the remainder of the paper, we will claim that 
"leptogenesis predicts an observable BR{£j ^i7)", if the Vl elements required exceeed the last column 
of table l2.lU last three rows), with y^ — 1- It is clear that the branching ratios can be decreased, for fixed 
Vl, by decreasing y^ and adjusting weak scale SUSY parameters. However, if SUSY is discovered, these 
masses and mixing angles could in principle be measured at colliders, and some information about the 
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magnitude of 2/3 could be available through the renormalisation group equations |45j . This assumption 
of universal soft masses at the scale Mx will not be crucial for our conclusions. Additional contributions 
to the off-diagonal soft masses are unlikely to cancel the ones we discuss, so the lower bounds we set on 
LFV branching ratios, from requiring leptogenesis to work, should remain. 



BR{^i ej) < 1.2 X 10"" 
(PSI) 


y3^L*32^L31 < .006 


BR{t ^ 67) < 2.7 X 10-'' 
(CLEO) 


y3^L*33^L31 ^ 12 


BR{t At7) < 1.1 X 10"'^ 
(CLEO) 


2/3^L*33^L32 < 9 


BR{fi -> 67) - 10-^*^^^ 
(PSI/nufact) 


y3^L*32^L31 - 3 X 10-^ 


BR{t 67) - 10"^ 
(BABAR/BELLE) 


y3^L*33^L31 ^ 0.2 


BR{t ^7) ~ 10-y 
(BABAR/BELLE/LHC) 


yiVE33VL32 - 0.2 



Table 1: Current limits 0^1 and hoped for sensitivities 02] of some experiments. The numerical bounds 
in the right column are multipHed by (^f--)^ (^) . If tan /3 is large, these rare decays are sensitive 
to smaller angles in Vl EH- 

Various CP violating phases in the neutrino and slepton mass matrices could be measured in 
upcoming experiments. However, the experimental sensitivity to the phases depends on the magnitude 
of unmeasured real parameters. Anticipated QvP(3 experiments may be sensitive to a maximal phase </>', 
for the neutrino mass spectrum we consider. The minimum value of the angle 8 that could be measured 
at a 1/ factory depends on Am32,Am|;^ and 6'i3 (see e.g. 021); so there is no forseeable clear upper 
bound. The imaginary part of the product of off-diagonal slepton masses 3{[TO?]i2[m?]23[TOp]3i} = 
J{'rfi\ — m\){m^ — in^){'m\ — ?7i|) could be measured in slepton flavour oscillations down to J = 10~^ 
|49j . J depends on the magnitude of the [m1]ij as well as their phases, so all the phases in the (s)lepton 
sector can be of order 1. 

2.2 Review of the parametrisation 

It was shown in |10| that the seesaw can be parametrised from the bottom-up, using the neutrino and 
sneutrino mass matrices. See |5()| for applications. A similar phcnomcnological parametrisation of the 
non-SUSY seesaw could be used, if the scale M of right-handed masses was low enough to measure 
dimension six operators oc 1/M^. We briefly review jJOl here. 
It is in principle possible to extract the matrix 

p = yIy,^vIdI.Vl (11) 

from its contribution to the renormalisation group running of the slepton mass matrix. This relies 
critically on having universal soft masses at the GUT scale, and on very precise measurements of 
sneutrino masses and decays. It is therefore unrealistic ^^l- However, since SUSY has not yet been 
discovered, Dy and Vl can be used as inputs in a "bottom-up" parametrisation of the seesaw. 
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The aim is to determine and M from [rrii,] and P. Vl and Dy can be determined from P, and 
used to strip the Yukawas off [rrii,]: 

D-Wl^^VlD-' = V^D-^V^ = M~\ (12) 

where the left hand side of this equation is known ([m^] is one of the inputs, and and Dy were 
obtained from eq. Hll|) '). Therefore, Vr and Dj^ can also be determined. This shows that, working 
in the basis where the charged lepton Yukawa coupling, Yg, the right-handed Majorana mass matrix, 
M, and the gauge interactions are all diagonal, it is possible to determine uniquely the heavy Majorana 
mass matrix, Ai, and the neutrino Yukawa coupling, Y^ = V^DyVL, starting from [m,^] and YjYi,. 



2.3 Leptogenesis, and the upper bound 

The see-saw mechanism provides a natural framework to generate the baryon asymmetry of the Universe, 
defined as rys = {ns —ng)/s, where s is the entropy density. As was shown by Sakharov|54j. generating 
a baryon asymmetry requires baryon number violation, C and CP violation, and a deviation from 
thermal equilibrium. These three conditions are fulfilled in the out-of-equilibrium decay of the right- 
handed neutrinos and sneutrinos in the early Universe. In the remainder of this paper, "right-handed 
neutrinos", and the shorthand notation i^n, refer to both right-handed neutrinos and right-handed 
sneutrinos. 

In gravity mediated SUSY breaking scenarios, these is an upper bound from gravitino production 
on the reheat temperature Tfljf of the Universe after inflation. The gravitino has a mass m.3/2 ~ "nisusY 
and only gravitational interactions with SM particles, so it is very weakly coupled, and long-lived. If a 
significant number of them decay at or after Big Bang Nucleosynthesis, they could disrupt the predicted 
abundances of light elements [Z|. Gravitinos can be created by various mechanisms in the early 
Universe, such as scattering in the thermal plasma^ d j or direct coupling to the inflaton (preheating) 
|55[ I56j. The latter is effective, but avoidable j56| . The number density of gravitinos produced in 
scattering increases with the plasma temperature, so the bound on 713/2 sets an upper bound on the 
reheat temperature of the Universe after inflation of 

Trh ^ 10^ -> lO^^GeV (13) 

(corresponding to 7713/2 ~ 100 GeV — > 10 TeV 7 ). This bound assumes that the gravitino decays; there 
are models where the gravitino is the LSP, which allow Thh ~ 10^^ GeV j57j . 

Let us briefly review the mechanism of generation of the B AU through leptogenesis ^ ■ At 
the end of inflation, a certain number density of right-handed neutrinos, n^^, is somehow produced. If 
these right-handed neutrinos vr. decay out of equilibrium, a lepton asymmetry can be created. The 
subsequent ratio of the lepton excess to the entropy density s is given by 

VL = = > ^ (14) 

S — ' s 

i 

The CP-violating parameter is determined by the particle physics model that gives the masses and 
couplings of the vn. The value of Uuj^/s depends on the mechanism to generate the right-handed 
neutrinos. We assume the vn-i are generated by Yukawa scattering in the thermal plasma, in which 
case Tii/^ /s ^ n'^'^/s ~ -2/5*, where n^'^ is the equilibrium number density of massless particles, and 
g^, ~ 230 is the number of propagating states in the supersymmetric plasma ^. This also implies an 

^in our conventions, n^^ = (n^^ + nj7^)/2. The .2 is an approximation to g^n'^'^ /s = f (3)135/(871*) . 
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upper bound on the mass: Mi ^ Trh- Finally, di is the fraction of the produced asymmetry that 
survives after Vfi decay. To ensure c?i ~ 1, lepton number violating interactions (decays, inverse decays 
and scatterings) must be out of equilibrium when the right-handed neutrinos decay. In the case of the 
lightest right-handed neutrino vr-^ , this corresponds approximately to 

K^-^J^<1 (15) 

where H is the Hubble parameter at the temperature T, and Td^ the vr^ decay rate. There are 
two competing requirements on the vr parameters — the couplings must be large enough to produce a 
thermal distribution, but small enough that the vr decay out of equilibrium. Thermal leptogenesis has 
been carefully studied in JHI'*- The numerical results of suggest that n^^di/s < n|;^/s: either 

n^^ does not attain its equilibrium number density, or lepton number violating interactions wash out a 
significant fraction of the asymmetry as it is produced. Defining an effective light neutrino "mass" 

Ml Afi ^ ' 

nvi^di/s 10^'' is realised for |^ 5 x 10^^ eV ^ fhi < 10^^ eV. The precise numerical bound on rhi 
depends on Mi, and can be found in [TS| . 

For mi > 10^^ eV and Mi ^ 10^ GeV, the dilution factor di can be approximated |58[ I59| 

(17) 

s 6(7, + 1 

with K ~ 910mi/ey from eqn H15|l . This is a slight modification of the approximation, to ensure 
that it falls between the Mi = 10® GeV and 10^° GeV fines of HI], in the relevant range .001 eV 
^ mi ^ .1 eV. The exact numerical factor is important, because it is difficult to get a large enough 
asymmetry. Multiplying di by a factor of a few significantly increases the parameter space where thermal 
leptogenesis can work. The approximation (|17|l neglects the decrease in di at mi ^ 10"'' eV, which 
is due to underproduction of vr^ in scattering. This is reasonable, because mi > mi 0, and we take 
mi = m2/10. 

The last step is the transformation of the lepton asymmetry into a baryon asymmetry by non- 
perturbative B-|-L violating (sphaleron) processes (1!, giving 

= C77B_L = (3-9) X 10-", (18) 

where C = 8/23 in the Minimal Supcrsymmctric Standard Model. Big Bang Nucleosynthesis constrains 
rjE to lie in the range of eqn (|18|l . In a flat Universe, the CMB determines rjB — (0.75 — 1.0) x 10"^" 
|53j. The wider BBN range is used in this paper, because it is difficult to generate a large enough rjB- 
The CP asymmetry can be approximated as 



ei 



3 Ml 



Stt [Y,Y,^]ii 



^El-{[Y.Yt]?,}f§) (19) 
Im^[Y.i^Y/]iiL (20) 



*See l'2Ul for a detailed analysis of thermal leptogenesis at higher temperatures, including the effects of Uii2 and 
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if the lepton asymmetry is generated in the decay of the hghtest right-handed neutrino, and if the masses 
of the right-handed neutrinos are hierarchical ^. Were the asymmetry produced in the decay of i/fj^ or 
un^, it would depend on a different combination of couplings. 
It is straightforward to show |S] that if ei is written 

^^MimsSHMY (21) 

then eqn H20(l implies the upper bound Shmy < 1- The numerical results of |21[ I60| agree with this 
constraint. Using eqns H14(l and (|18|l . this can be transformed into a lower bound on Mi: 

Setting 7713 to its 90% C.L. upper bound 0.063 eV, and di to its maximum value n'^'^/s ~ A5/(2Tr'^g^), 
implies Mi > 3 x 10^ GeV. 

This lower bound on Afi comes very close to the gravitino bound eqn H13(l on the reheat temperature. 
Thermal production of the Vfi requires Mi ^ Trh so either e is close to its upper bound, or Mi, Trh > 
10^ GeV, or thermal leptogenesis does not generate the observed baryon asymmetry. We explore the 
first option, and somewhat the second. The third possibility, non-thermal i'r production, has been 
discussed by many authors (see e.g. references of 



3 1713 ~ 

a 2" ^1 



3 Analytic approximations for Shmyi Mi, ffii 

At least three inputs are required to parametrise thermal leptogenesis dHl im EHI- A possible choice 
would be the mass Mi and decay rate cx mi of the vr, and the CP asymmetry ei. However, e oc Mi, so 
we use Ml, rhi and Shmy (introduced by Hamaguchi, Murayama and Yanagida), where Shmy measures 
how close e comes to saturating its upper bound. Note however, that Shmy is not a CP phase. 

This section contains simple analytic approximations indicating the dependence of leptogenesis 
parameters on measurable quantities, such as neutrino masses and rare LFV decays. We used this ap- 
proximation, with attention to the phases, in '37!- A similar, somewhat simplified version was introduced 
in [ini. 

The inputs for the analytic approximation are: 

Vl, Dy, U, [m,] (23) 

Two of the angles of U are known, and the CHOOZ angle is bounded above. The eigenvalues y.i of the 
neutrino Yukawa matrix are unknown, and realistically cannot be determined from the sneutrino mass 
matrix. It seems reasonable to assume a hierarchy for the {yi}, since we measure hierarchical Yukawas 
for the quarks and charged leptons. The yi remain as variables in the equations; we will discover that 
only the smallest eigenvalue yi is relevant, and can be "traded" for the mass Mi of the lightest i/n, 
which is tightly constrained. Vl contains three unknown angles, related to the lepton flavour violating 
decays £j — > £ij. There are three phases in both U and Vl, all are unknown, and assumed to be chosen 
to maximise the baryon asymmetry. 

The lightest eigenvalue and corresponding eigenvector of M are estimated in the first Appendix, 
which also contains some simple (but illuminating) 3-d plots of leptogenesis parameters. The mass of 

^If the hierarchy in is similar to that of the quarks and charged leptons, then a hierarchy in the Mi is natural. 
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the hghtest vr is 

where the matrix W — VlU is the rotation from the basis where the vl masses are diagonal to the 
basis where the neutrino Yukawa matrix YJY^ is diagonal. There are three limiting values for Mi, 
corresponding to Mi ~ ylvl^/nn: Mi yfv^/m^ when W13 1, Mi ylvl^/mi when W13, W12 — > 0, 
and Ml y'lv'^/m^ when Wiz < m^/m^, W12 1. This is easy to see in figure^ 

Ml is the only quantity relevant for leptogenesis which depends on yi. The latter is effectively 
unmeasurable; it is constrained by theoretical expectations, and by the requirement that the analytic 
approximation be self-consistent. Theoretically, the eigenvalues of YJy are expected to be hierarchical, and 
of order the quark or lepton Yukawas, so figure^lis plotted with yi ~ 10~^. The approximations of this 
section are consistent, provided that the dropped 0{yf,y1/y2) terms are smaller than the O^mi/m^) 
terms which are kept. This is the case for yi ~ 10"'^. Since yi is unmeasurable and only weakly 
constrained, it can be adjusted, as function olrrii and Wij, to obtain a value for Mi where leptogenesis 
could work. In fact, since Mi cx yi is tightly constrained, the requirement Mi ~ 10^ GeV "determines" 

yi- 

The eigenvector (130)1 can be used to evaluate the uri decay rate: eqn. (fTC|l becomes 

fhi has various limits: rhi ms for W13 large, mi 7712 for Vl^i2 large and W13 < 7772/7773 , and 
7711 nil when W ^ 1. This is easy to see from the RHS of figure^l In the W ^ 1 limit, washout is 
minimised, because the dilution factor di oc 1/7/11. 

To saturate the upper bound H21|l . 6hmy needs to approach 1. Evaluating eq. (|20|l with the 
eigenvector (jSnj, gives 



\WiiWi2\^mim\ + \WiiWiz\^miml + \Wi2Wij,\^m2m\ 



(26) 



This paper is about the relation between real low energy observables (such as BR{fi — > 67)) and the 
baryon asymmetry, so scant attention will be paid to the phases in U and Vl- For most of parameter 
space ^, the phases can be chosen such that Shmy is larger than the second expression in eqn (|26|l . 

This second expression is plotted on the LHS in figure [3 Shmy can approach 1 if the numerator 
(upstairs) is dominated by 77737772 or by 77737771. This is because of the 7713 in the denominator. If the 77i|mi 
dominates upstairs, then Shmy will approach 1 when W^jirij ~ Wiimi and \Win\'^m'^ ~ |M^i3pw§, or 
equivalently, when 

< Wio < — and W12 < — , — I (27) 

mj 77I3 ^"^7712 7771 ^ ^ 

This corresponds to the highest ridge in Shmy in figure Notice that the position of the peak depends 
sensitively on the lightest neutrino mass 7771. 

^everywhere but when the three terms upstairs have equal magnitude 



9 



If 171^1712 dominates upstairs, then Shmy ^ 1 when 



t2 



W?^'^ < Wl, < W^2— and Wl, < W^^ — . W?,^ (28) 
TO3 7713 n^i n^i 

This corresponds to the shoulder at shghtly large W13, which is cut by 1^x2 ~ 1 in the LH plot of figure 

u 

Finally, for W13 very small, S 7^2/7)7,3 ^ 0.1 along the ridge at W12 ~ 0.1. This corresponds to 
the 7)7,2'77i term dominating upstairs, and arises when 

W?,^ < < and Wl, < W^,^, W^^^ (29) 

7772 "^2 ni3 7773 

Although Shmy does not reach its maximum value for these parameters, the washout is small, so the 
baryon asymmetry generated is only slightly too small. As we will see in figure |21 it is large enough if 
Mi.Trh ~ 10^° GeV are allowed. 



4 When does thermal leptogenesis work? 

The baryon asymmetry can be approximated as 

8di 3 , 
■HB ^ -jrir^ — jMim,3dHMY (30) 

by combining eqns lfn|) . ((TH|) . and (|^ . This is plotted in figure Q] which suggests that t/b can be large 
enough. 

The issue is whether a large enough asymmetry can be generated, so the observational upper 
limit on tjb is unimportant. Also, the asymmetry calculated here is the upper bound corresponding to 
maximal CP violation, so it can be reduced by taking smaller phases. We use the one-cr observational 
lower bound on rjB from nucleosynthesis: 7/3 3 x 10~^^. To obtain a large enough baryon asymmetry 
by thermal leptogenesis, the parameters Mi, Shmy, and di must occupy narrow ranges. The washout 
effects are minimised when the vj^ decay rate is small, which corresponds to 14^ 1. In this case, 
n-undi/s = di < 10~^ which implies the lower bound e lO"*". (If e <; 10~^ can be obtained, then 
di ~ 10~^ is large enough.) Eqn (|21|l implies a lower bound on Mi to get ei large enough. In addition. 
Ml ^ Trh is required for thermal production; the canonical SUSY gravitino bound is Trh ^ 10^ GeV, 
so 5e X 10^'''GeV < Mi < Trh- Since e ~ 10"'^ is required, for Mi <, 10^ GeV one must have Shmy 1- 
The parameter space of choice can be summarised as 

few X lO^Gey 
di 

Shmy 

As can be seen from figures 01 — El it is difficult to simultaneously satisfy these conditions. Mi and di 
increase as W13, W12 0, but S decreases. 

The analytic approximations of the previous section show that the baryon asymmetry depends on 
U and the first row of Vl (via Wij), on the light neutrino masses 771^, on the lightest neutrino Yukawa 
7/1, and on phases. These real parameters are known, or could be experimentally constrained in the next 



£ Ml < few X lO^GeV 
45 

-> 1 



(31) 
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Figure 1: 3-d plot of tje = 8die/23, for central neutrino mass values, toi — m2/10 and Mi = 10^ GeV. 
On the left, rjB is plotted as a function of u)i2 — log[Wi2] and W13 ~ Zo5[iyi3]. On the right, ?]b is 
plotted as a function of ^13 and xi2, defined such that Wu = cos6\Yi3sin{9soi — 10'^i27r/2). The RHS 
measure on parameter space is more sensible, see the discussion after ean (|34|l . 

20 years — with the exception of mi, j/i, and Vli2- For the purposes of this paper, Vli2 is included with 
the measurable angles, and yi is determined as a function of mi, by requiring that Mi be in the range 
H31|l where leptogenesis could work. The baryon asymmetry then becomes independent of yi. Some 
subtle dependence on mi remains: the area and location of the high ridge in figure^ depend on mi, but 
the baryon asymmetry and low energy footprints do not. This is discussed in the Appendix about mi. 

Notice also that it could be expected to have a similar hierarchy in the neutrino Yukawas as in the 
other fermions, in which case yi ~ /i^, hg or h^- This gives 

where W^^mj is usually of order m2- If yi — h^, then 773 is too small over most of parameter space. 
This was found in some models by jl5| . The baryon asymmetry can be large enough, for yi ~ /i„, in 
the small area of parameter space where Wi^mj ~ mi ^ m2/100. This is in the mi Appendix too. 

The baryon asymmetry depends weakly on tan/3, when Mi is taken as an input, and di is ap- 
proximated as (X 1/mi cx sin^ /?. The m^ are experimentally measured, and therefore independent 
of sin^ so it is clear from eqn (|3()|l that the sin/3 dependence arises entirely from mi. If instead 
Ml = (jjfv'^)/\Wijmj\, then rjB cx sin^/3. In both cases, larger sin/3 is marginally favoured. 

The parameters W12 and W13 are convenient, because they summarise the unknown mixing angles 
and phases. The physically relevant quantities for leptogenesis (Mi, e,...) can be plotted as a function of 
the two real unknowns \ Wi2\ and |W^i3|. However, the Wij are not observable — the matrix W is related 
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to the more physical matrices Vl and U hy W = VlU. Recall that Vl rotates from the basis where the 
neutrino Yukawa matrix is diagonal to the basis where is diagonal, and U rotates from the basis 
where [m^] is diagonal to the basis where Ye is diagonal. W can be written 

Wi3 = VLiisineise-'^ + VL12/V2 + VL12/V2 (33) 
1^12 = Vlii sindsoi + Vli2{cos 6 sol - s'md soi sin 6*136"^)/ V2 

- 13 (cos Osoi + sin Osoi sin 013 e'-^ ) / 72 (34) 

where 9i2 — Ogoi and 623 = 7r/4 in the MNS matrix. 

From a model building perspective [3T], there are two natural limits for W. The most popular is for 
the large leptonic mixing angles to come from the seesaw structure of the light neutrino mass matrix. In 
this case, Vl ^ I can easily arise, soW^U. This is similar to the quark sector, where the CKM matrix 
(the analogue of Vl) has small angles. An example of this is texture models where the large atmospheric 
mixing angle is due to a Vfi mass eigenstate having approximately equal Yukawa couplings to and 
Vr |S2] ■ Alternatively, the electron Yukawa matrix Ye could be "the odd man out" ; it could have large 
off-diagonal elements in a basis where the neutrino mass matrix [m^] and YJy are simultaneously almost 
diagonal In this case Vl ~ and W ^ 1. These two cases are discussed in the following two 
subsections. Figure ^ suggests that thermal leptogenesis can work for W ~ 1. As we shall see, a large 
enough asymmetry is also possible in the Vl — > 1 limit, if a slightly larger T^h is allowed. 

The plots are functions of log W13 and log W12, rather than, e.g.Wi2 and W13. It is sensible to use 
logarithmic measure on unknown physical parameters^ because it is equally probable for mixing angles 
to have any order of magnitude between e.g.lO~^ and 1. However, Wij are not physical parameters, 
so this reasoning does not apply to them. Specifically, values of W12 ^ sinOgoi arise in the presumably 
small area of parameter space where Vl — U'^ . This reasoning does apply to the CHOOZ angle and the 
unknown angles of Vl, but we prefer to plot rye as a function of two unknowns, rather than four. So 
a more appropriate measure on W12 might be logarithmic in the difference away from U12 — sin Oggi . 
Therefore, on the RHS of figuren]is plotted the same function as on the LHS, but as a function of Wis ~ 
log(Wi3), and xi2, the latter defined such that 

1^12 =sin(0,„,- 10X1^^2). (35) 

X12 ~ log[Wi2 — Ui2\ is an approximation to the log of the unknown angles of Vl (see eqn (I34II 1 . So the 
RHS plot tells us the same information as its twin on the left: the asymmetry is largest if a large angle 
in Vl cancels the large solar angle in the MNS matrix U . 

A final technical comment: W ^ 1 and Vl ^ ^ mean the 12 and 13 matrix elements are small 
<^ .1. W ^ 1 means that W maximises 773, so W12 ^ .1 and .01 ^ W13 ^ .1. Vl ^ I includes both the 
possibilities that the angles of Vl are smaller, or larger, than the CHOOZ angle. 

Section UTI studies the phenomenological consequences of sitting in the region where thermal lep- 
togenesis works easily, which corresponds approximately to W12 ^ -l, .01 ^ W13 ^ .1. Then in section 
14.21 some of the parameters which are fixed in figure ^ are varied, so a large enough baryon asymmetry 
can be generated for Vl — I- The parameter space between these two limits is discussed in section IHH 

4.1 Wr^l 

The parameter space where 775 is largest in figure ^ corresponds roughly to 

.01 ^ Wis ^ -1 W12 ^ .1 . (36) 

'^this choice of measure is neither unique nor universally agreed on 
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This can be understood from the analytic approximation (j^OJ). We fix Mi ~ 10^ GeV, so r]B oc (IiShmy- 
The factor di is largest when the vr decay rate T cx rhi is smallest, so more of the asymmetry survives 
when W ^ 1 (see the expression 1^ . Shmy parametrises how close e can come to its upper bound 
iPT|l . For nil ~ TO2/IO, the values of W^i2, W13 where 6hmy is maximised fean l77|) correspond to eqn 
H35|l . riB is maximal at smaller W13 than 6hmy, as can be seen by comparing figures and 01 This is 
due to the lepton number washout encoded in di, which is faster at larger W13. 
To obtain W12 and W13 in this region, Vl must have the form 

Vl=R23[U^]^ (37) 

where R23 is an unspecified complex rotation in the 23 plane (written in the form of eqn © with 
012 — ^13 — 0, S = sin 6*23, C = cos 6*23, and taking a 23 phase a), and [U~] is a matrix whose angles are 
roughly those of the MNS matrix ±.1. The unknown i?23 appears because leptogenesis only depends on 
the first row of W. 

It is interesting to study the implications for £j £ij of eqn (|37|l . Taking [U~] = U 
Vl31 = 5e*("+^'/2)ci3Si2 + Csi3e'* 

VlZ2 = 5e^("+^'/2) (C23C12 - S23Sl3Sl2e'') + Cs23Cl3 

^ 5e^("+^'/2)c^^^/V2 + C/V2 

For generic values of 5, this implies VL-i2 ~ 1, so an experimentally accessible t ^ ji'y branching ratio. 
If, on the other hand, S is tuned to make Vl32 — > 0, then BR(t — > //7) would be unobservable. However, 
in this case Vlzi ~ — sin 9 sol + cos^ 9 sol — 1/V3, so r ^ 67 should be observable. One can conclude 
that if leptogenesis takes place in the W ^ \ peak of figure then one or both oi t ^ fi'y and r 57 
should have a branching ratio ^ 10~^ (according to the leading log approximation of the introduction). 
Similarly, BR(pL 67) should be > 10-^'* if S or S13 Z IQ-^. 

4.2 Vl = 1 

It is barely possible to get a large enough baryon asymmetry when the angles in Vl are small, although 
this is not evident from figure ^ In the limit Vl 1, the matrix W ^ U, so W12 — sin 61^0; and 
Wi3 ~ sin 013. In figure 775 is at least a factor of 6 to small at log 1^12 ~ —0.5 (equivalently, X12 
small), but there is a bump at W13 '--^ .1. In this section, rjB at Vl = 1 is increased by varying m^, m2 
and Afi; Vl close to the identity is discussed in the following subsection. 
If Vl = 1, eqns and ((^ give 

M, = J^- 

TO2SJ2 

rhi — TO2 

To maximise the asymmetry. Mi is taken to be / x 10^ GeV, where / is a few (this determines j/i — 
Mim2s\2/v'^ ~ 7.2/ X 10~* ). In figure 10 of JT], the lepton asymmetry is plotted as a function of fhi 
for various values of Mi and ei = 10~^. This plot shows that for e ~ 10~^, a large enough asymmetry 
can be generated when fhi ~ m2. Note that the washout effects are correctly included in this plot of 
|11|. so this result does not depend on the analytic approximation of eqn H17() . Also, the asymmetry 
increases as fhi decreases, so small values of the solar mass are prefered. 
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The upper bound of eqn implies that for e = 10~^, fm^SHMY = 0.25eV. From the experimen- 
tally allowed range of TO3 given after eqn we see that A/i ~ 3 x 10^ GeV is required, assuming 5 ~ 1 
is also possible. For Vl — I, 

sfamj^gf^mz sinjcj)' - 25) 
Ohmy — r-2 —. — 2 2\"^ 

which approaches 1 when S13 ~ S12TO2/TO3. In the RH plot of figure (Q), Vl = 1 corresponds to 
Wis, — sin 013 and X12 —00. In figure|21 &hmy is plotted as a function of ^13 on the LHS; the analytic 
approximation to r]B (eqn H3U|l 'l is plotted on the RHS. So thermal leptogenesis "works" at Vl = 1, for 
Ml ~ 6 X lO'^ GeV. 

Phrased another way: for arbitrarily small ij — > ^^7 branching ratios, requiring thermal leptogenesis 
to work predicts the CHOOZ angle 6*13. If 7713 is taken at its 90%C.L. upper bound, and Mi ~ 6 x 10^ 
GeV, then rjB ^ 3^ 10~^^ can be obtained. This requires a CHOOZ angle of 6*13 '--^ 4 x 10~^, and phases 
which satisfy 2S — (p' — Tr/2. 




2e-1 1 



-1.5e-11 



1e-11 



-5e-12 



Figure 2: On the LHS 6hmy, and on the RHS the analytic approximation to rjs, as a function of 
513 = log[sm9i3]. Shmy ^ 1 is required to get a large enough asymmetry, see the discussion in section 
14.21 The remaining parameters are tan^ 9soi — 0.44, TO2 — 7 x 10~^ eV, 7713 ~ 6.3 x 10~^ eV, and for 
the r]B plot, = 1 and Ml = 4 X 10^ GeV. 



14 



4.3 Vl from 1 to f/t 



It is clear from the RHS of figure ^ that if a large enough asymmetry can be generated at = 1 
(Xi2 ~ —3), then enough baryons can be generated along the ridge leading to the peak, and also along 
the "other" ridge at W12 ~ mi/m2- To sit on this second ridge requires ~ 1, so has the same 
experimental signatures as discussed in section 14.11 This section is about the W13 ^ m2 /m^ ridge 
stretching from Vl ~ 1 to the W ^ 1 peak. 

Starting from the small X12, Hat section of the ridge and moving towards the peak corresponds to 
allowing small matrix elements Vli2, Vlis ^ .1- In this limit, 

Wi3 ^ sin6'i3 + FLi2/%/2 + Vli3/%/2 

and Wi3 ~ 0.04 is required to get a large enough asymmetry. Unfortunately, Wi^ ^ 0.04 determines 
a sum of three unknowns: ^13 ^ .04 could be observed, Vli3 ^ 0.04 induces a potentially observable 
r — > 67 signal, but Vli2 has no observable consequences. For Vl ~ 1, the Vli2 contribution to to? (eqn 
©) is suppressed by ?/| ~ 10^*. 

Figure|31is a contour plot in uis and X12 space of the approximation pOI) to rjB- The contours enclose 
the area when rjB > '^x 10~^^, for Mi = / x 10^ GeV, central values of TO3 and TO2, toi = TO2/IO, and are 
labelled by /. Allowing / > 1 significantly increases the available parameter space. This corresponds to 
increasing Mi (which should be ^ Trh), or increasing TO3, which is constrained by atmospheric neurtino 
oscillations, or for xi2 ~ 0.5, to decreasing TO2, which is constrained by solar neutrino experiments and 
KamLAND. Perhaps the most palatable way to increase / is to allow Thh ^ 10^° GeV. The value of 
riB chosen, rjB = 2 x 10""'^'^, is minimal. To obtain the CMB favoured 77^ ~ 9 x 10~^^, would require 
values of / that were four times larger. 

5 CP violation 

In a previous paper |37j, we discussed the relation between the leptonic phases that could be measured 
at low energy, and the CP violation required for leptogenesis. We assumed that e was large enough, 
and studied the relative importance of the neutrino factory phase S and the double beta decay phase (f)' 
for leptogenesis. If the right-handed neutrinos vm are produced non-thermally, getting e large enough 
may not be a significant constraint (see e.g. |BJ for a discussion and references). However, we have 
seen that it is a challenge when the vri are produced thermally. So in this section, we briefly discuss 
the relative importance of low-energy phases for thermal leptogenesis — imposing the constraint that e 
is large enough. 

It is well known that there is no linear relation between the "leptogenesis phase" and S or (j>' |64| . 
That is, the lepton asymmetry can be non-zero when S = (j)' = 0, and it can be zero when 5, 4>' 7^ 0. To 
overcome this, we introduced a statistical notion of "overlap" between the leptogenesis phase and the 
low energy phases of our parametrisation. The overlap Og aimed to quantify the relative importance of 
the phase 5 for leptogenesis, assuming that all the low-energy phases were 0(1). In [37], we considered 
the cases where Wi3Wi2'mlm2, or 14^2^2 ^11 '^2"^! ; is the most important term upstairs in 5hmy- That 
is, we consider Vl = I, Vl ^ I and the case of large Vl angles that do not exactly cancel those in U. 
This occurs over most of the parameter space where e could be large enough. However, e is largest in 
the small area of parameter space where W ^ 1 and Wi^Wiimlnii dominates upstairs in Shmy- So 
let us now consider which low-energy phases are important for leptogenesis in this case. 

Writing the phases explicitely gives 

e CX ^{WiiW^i} = ^{e'^[VLnCi3Ci2 + \VL12\e'^'^'K-C23Si2 - S23Cl2Sl3e*^) 
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Figure 3: Contour plot of tjb, as a function of oji^ ~ log[VFi3] and Xi2 — log[VLi2 + Vlis]. The contours 
enclose the area when r]B > >^ 10~^^, for Mi = / x 10^ GeV, central values of ms and TO2, and 
mi = 7712/10. In the direction of increasing area, the lines correpond to / = 1, 3, 6 and 9. 



+ |FLi3|e*2^-(s23Si2 - C23Ci2Si3e*^)]2 X [yLiiSise"*^ + |Fii2|e*2v-s23c,3 

+ |VLl3|e'2^^^C23Cl3]2} (39) 

where Lpij is the phase of Vij. Vij ^ ?7*ie*"i, because W ^ diag{e'^'^^ , e^'^^ ,1} ^. The "ncutrinolcss 
double beta decay phase" </>' is irrelevant for leptogenesis, because it only enters into M^i2- The phase 
of mi will always be important, because Wn oc e~"^/^, so e will be a sum of terms oc sin(m0 + ...). 
Both the phases (pi2 and (pi^ of Vli2 and Vli3 arc likely to have significant overlap with the leptogenesis 
phase, because Vl — so the \VLij\ are large. The "neutrino factory phase" 6 always multiplies the 
CHOOZ angle, which suppresses its contribution to e. 

The three weak-scale phases which have significant "overlap" with the leptogenesis phase, in the 
area of parameter space near ~ 1, are therefore (pi2, (fis and (fi. This is unfortunate, because although 

®This constraint on W is what we usually refer to as W ~ 1. Since Vii is real, uji ~ —0/2. 



16 



there is some hope of measuring and 5, there is no foreseeable experiment to determine any of these 
three. 

6 Discussion and summary 

This paper has discussed leptogenesis in a minimal model of the SUSY seesaw, with gravity mediated 
SUSY breaking and universal soft masses at a high scale. It uses a parametrisation of the model in 
terms of 

D,„^,U,Dy,Vl (40) 

where is the diagonal light majorana neutrino mass matrix (assumed hierarchical), U is the MNS 
matrix, Dy is the diagonal neutrino Yukawa matrix, and Vl diagonalises Y^Y^. The notation is briefly 
defined in table, giving the equation numbers of more detailed definitions. The angles of the unitary 
matrix Vl can be related, in SUSY models, to the rates for ij because Y^Y^ contributes to the 

renormalisation group equations for the slepton masses. This allows the right handed neutrino masses 
and Yukawa couplings to be expressed as a function of quantities which could be measured, in principle 
or in practise, at the weak scale. This parametrisation is briefly reviewed in section 



Y.^vIDyVl 


neutrino Yukawa, eigenvalues 


m 


M,Mi 


mass matrix, lightest eigenvalue 




Vl 


VlY,^Y,vI = Dl 




ipij 


phases of Vl 


sectional 




light neutrino mass matrix 


11 


U 


MNS matrix 


m 




angles; phases of U 


Em 


W 


VlU 


m 




^log [Wij] 


ED 


X12 


~log [U12 - W12] 


HI 


VL 


Icpton asymmetry 


m 


e 


CP asymmetry 




Shmy 


^/ ^max 


121126I 


fhi 


(X VRi decay rate 






dilution factor of lepton asymmetry 


1171141 


VB 


baryon asymmetry 





Table 2: Table of notation, with a brief description and the equation number of a more complete 
defintion. 

The baryon asymmetry produced in leptogenesis depends on the number density of right-handed 
neutrinos which decay, on e = the average lepton asymmetry produced per decay, and on the survival 
probability of the asymmetry in the thermal plasma after it is produced. We consider the "thermal 
leptogensis" scenario, in which the right-handed neutrinos are produced by scattering interactions in 
the plasma. Non-thermal production mechanisms are also possible, perhaps even probable, but depend 
on additional parameters from the sector which produces the right-handed neutrinos. Both the thermally 
produced Vfi number density, and the survival probability of the lepton asymmetry in the plasma after 
it is produced, can be computed from the reheat temperature of the plasma after inflation Trhi 
from the seesaw parameters. These processes have been carefully studied in |18llll| . 
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A convenient analytic approximation to the numerical results of JHI is used in this paper. A single 
function di (see equation ((T7|l ') is defined as the number density of i^i^x the survival probabilty of the 
lepton asymmetry once it is produced. So the baryon to entropy ratio today is r]B — 8die/23. ^ 

The CP asymmetry produced in the decay of the lightest is bounded above (for hierarchical Mi 
and TTij): 



where Mi is the mass of the vri and — •^/Arof^. Since di cannot exceed 45/(27r'*g*) for thermally 
produced Vfn, obtaining rys > 3 x 10~^^ requires Mi > 3 x 10^ GeV. 

In section|31 approximate analytic formulae for the lightest mass Mi , for the CP asymmetry e and 
for the baryon asymmetry r]B, are given in terms of our weak-scale parameters. These approximations 
are valid for hierarchical Mi and neutrino Yukawa eigenvalues yj . The baryon asymmetry can be written 
as a function 

?7_B(m2, r7i3, 6*23, 6*12; 6*13, Vli2, Vlis, yi, mi, phases) (42) 

where "known" low-energy parameters precede the semi-colon. We concentrate on the dependence of 
r]B on real parameters, assuming that the phases can be chosen to maximise the asymmetry. 

It is interesting that the baryon asymmetry only depends on 5 of the 8 unknown real parameters in 
eqn Two of these, 613 and Vli3, are possibly measurable; the constraints that thermal leptogenesis 
imposes on them will be discussed later for different areas of parameter space. On the other hand, there 
are no foreseen experiments that could determine yi, mi, and Vli2- Vli2 is included in the discussion 
with 6*13 and Vlis, because these three unknowns can be exchanged for the 12 and 23 elements of 
W = Vi,U . This simplifies expressions and is convenient for plotting. The dependence of rjB on mi is 
subtle, comparatively unimportant, and discussed in an Appendix. The lightest right-handed neutrino 
mass Ml, and therefore the baryon asymmetry, is proportional to y\. So y\ is exchanged for Mi. 
This is a peculiar exchange — why do we want to use a GUT-scale mass as input in our weak-scale 
parametrisation? The off-diagonal elements of Vl, and the yi, are related to lepton flavour violating 
off-diagonal slepton mass matrix entries (to which processes like ij iij are sensitive), and to slepton 
mass differences. The smallest neutrino Yukawa yi makes negligeable contributions to both these effects. 
However, Mi > 3x 10® GeV is required for thermal leptogenesis to have any hope of working, and if SUSY 
is discovered, the sparticle spectrum could give some indication of the gravitino mass, and therefore the 
allowed reheat temperature Trh > Mi. So we "determine" yi by requiring that thermal leptogenesis 
could produce a large enough asymmetry: 3 x 10® GeV < Mi < Trh- Then we study the requirements 
on the remaining parameters such that the asymmetry is large enough. These additional requirements 
may have observable consequences. 

The analytic formulae of section |31 are simple and compact, but nonetheless difficult to visualise. 
The asymmetry depends on the first row of the matrix W , so for qualitative understanding, we show 
3-dimensional figures of leptogenesis parameters as a function of 0^13 ~ \0gW13, and UJ12 — logVFi2. 
Logarithmic measure is reasonable for unmeasured but observable matrix elements — which the Wij are 
not. For small angles in Vl (see section 0] for a general discussion), the Wij can be related to the more 
physical matrix elements Vliu and Uij-. W13 ~ 9i3 + Vli2 + Vli3, W12 sin6'so; + VLi2 + VLi3- To present 
the area of parameter space where leptogenesis works with a more physical measure, we therefore plot 
the baryon asymmetry as a function of UJ13 and X12 ~ log[VLi2 + VL13] in figures ^ and |31 

We define thermal leptogenesis to "work" if it can produce 773 ^ 3 x 10~^^, as required by Big 
Bang Nucleosynthesis. For Mi ~ 10^ GeV (consistent with the canonical gravitino bound T^b ^ 10^ 

9 where the 8/23 arises in the transformation of the lepton asymmetry into a baryon asymmetry by the electroweak 
B + L violating processes. 
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GeV), leptogenesis can work: there is a limited parameter space where the upper bound on e is almost 
saturated, and di is close to maximal. This can be seen in figure 13 where the baryon asymmetry is 
large enough inside the contours, which are labelled by /, where Mi = / x 10^ GeV. Increasing Trh 
(and thereby the allowed Mi) enlarges the parameter space where thermal leptogenesis works. 

We now come to the aim of the paper — what are the weak scale foot prints of thermal leptogenesis? 
What parameter values must be observed, if thermal leptogenesis works in an MSUGRA model? 

Suppose first that / ~ 1, which corresponds to Mi ~ 10^ GeV for central values of the light neutrino 
masses. Thermal leptogenesis works in the area of figure Oat ojis ~ —2 and Xi2 ^ —0.5. This small 
area of parameter space is discussed in section lOl and occurs if = VlJ7 ~ 1. The phenomenological 
consequences of this area of parameter space are unambiguous: the branching ratio of r fij, or 
r 67, should be large. More concretely, at least one of Vl32 or Vlsi is 0(l/\/2), so according to the 
estimates of table O BR{t £7) > 10"^ 

From a theoretical model building perspective, this area of parameter space corresponds to the 
neutrino Yukawa and light mass matrices Y^, and [rriu] being almost simultaneously diagonalisable. The 
large leptonic mixing angles arise in the rotation from this basis to the one where the charged lepton 
Yukawa matrix Yg is diagonal. 

The baryon asymmetry is largest at this point for two reasons. The vn decay rate (eqn is 
small, so lepton number violation is slow after the asymmetry is produced, and more of the asymmetry 
survives. Secondly, the asymmetry produced is almost maximal; it comes within a factor of 0(1) of the 
upper bound eqn (|^ . This is discussed after eqn 

The area of this parameter space, where rjB is largest, depends on the smallest neutrino mass mi. 
The plots are made with mi = m2/10; the area shrinks as mi decreases. This peak in rjB only exists for 
mi ^ 0. It is interesting that the CP asymmetry and the low-energy footprints of this area of parameter 
space are independent of mi. However, the number density of ua (and therefore the baryon asymmetry) 
decreases for mi ^ 10~^ eV, and our approximation fails. See the Appendix for a discussion. 

In brief, if a sparticle spectrum consistent with gravity mediated SUSY breaking is measured, with 
a gravitino mass of m3/2 ~ 100 GeV (Trh ^ 10^ GeV), then r ^7 or r ^ 67 must be observable for 
thermal leptogenesis to work. 

Now consider the enlarged parameter space allowed by Mi/ (10^ GeV) = / > 1 in figure^ thermal 
leptogenesis works for W13 ^ 1112/ 111,3 and pretty much all values of 14^12. This sets one constraint on the 
three "physical" matrix elements sin 6*13, Vli3 and Vli2- As discussed in sections l4 . 21 and f4 . 31 it can be 
satisfied if any one of the angles is 0(m2/m3). These possibilities have different weak-scale implications. 

If Vl has small angles like the CKM matrix, Vli3,Vli2 ^0.1, then leptogenesis requires that the 
CHOOZ angle ^13 ~ 0.04, which is close to its current experimental bound. This implies that the baryon 
asymmetry is determined by parameters which can be measured in the neutrino sector m3, m2, ^13 
and the phases d and cj)' . 

If the CHOOZ angle ^13 <^ 0.1, then it is still possible to sit on the W13 ~ 0.04 ridge, by having 
yLi3, or Vli2 ~ 0.04. The former angle is related to r ^ 67, and could perhaps be measured in this 
process. Unfortunately, Vli2 appears in the RGBs multiplying y^, the middle Yukawa eigenvalue, so 
has no observable consequences in the slepton mass matrix. So thermal leptogenesis can "work" when 
^13 and BR{lj £^7) are unobservably small. 

The matrix elements Vlisi and Vli2 are small along most of the W13 ridge currently under dis- 
cussion. Many texture models occupy this area of parameters space, where the CKM-like matrix Vl 
(between the bases where Y^, and Ye are diagonal) has small angles. The large angles of the MNS matrix 
then arise from the majorana structure of Y^ M^^Y . 

^"Caveat: rjg oc Mim-^ / iri2 in this case, and we "determine" Mi ~ 6 X 10^ GeV by requiring rjg large enough. If m,3 
(m2) is larger (smaller) than the current best-fit values, 77s increases. 



19 



The baryon asymmetry is larger along the ridge than in the rest of parameter space, because the 
washout is moderate — fh ~ m2 — and because the CP asymmetry e approaches its upper bound 1)4111 . 
Notice that the baryon asymmetry on this ridge is independent of mi — the solution remains as mi 0. 
As discussed after eqn (|26f) . there are two limits in which e is maximal: the ridge where W13 ~ 7712/7713, 
and the previously discussed peak. There is an orthogonal ridge in figureOl at xi2 ~ —0.5 {W12 ~ 0.1), 
where Shmy ~ '712/^13, but washout is minimised. It has the same observable footprints as the peak. 

Until now we have considered if thermal leptogenesis works, then what should we see at low en- 
ergy? Allowing TfiH ^ 10^° GeV, it seems just about all phenomenology is consistent with thermal 
leptogenesis: large r jij, 6*13 ^ 0.04, observable t ej, nothing observable at all. ..So now consider 
the inverse question: are there weak-scale observations that can rule out thermal leptogenesis? The 
previous discussion is vague because SUSY has not been discovered. Clearly thermal leptogenesis does 
not work if W13 is too big or too small. Since all the terms which contribute to W13 cannot be measured, 
no experimental lower bound can be set. However, one could tell that W13 is too large, for instance if 
large Vli3 (t — > 67) is measured 

The analysis of this paper relies crucially on the assumption that the vr are produced thermally. 
A larger number density of the lightest vr, n^j^/s, could be produced non-thermally, so a large enough 
baryon asymmetry could be produced with a smaller e. This would enlarge the available parameter 
space. Furthermore, if the vr are produced non-thermally, they could be vr2 or vrz, making the 
formulae for e inapplicable. 

In summary, we study the baryon asymmetry resulting from the decay of the lightest right-handed 
neutrino Ur^ , assuming the vr-^^ s are produced thermally. We present compact analytic approximations 
for the quantities relevant to thermal leptogenesis, in terms of the light neutrino masses, the MNS matrix, 
the smallest eigenvalue of the neutrino Yukawa matrix YJy, and the matrix Vl which diagonalises Yi, on 
its SU(2) doublet indices. In the MSUGRA scenario, we can trade these parameters for the neutrino 
and sneutrino mass matrices (m^ and mf), or more usefully, for rui,, for the branching ratios of lepton 
flavour violating decays £j — > £ij, and for the lightest right-handed neutrino mass Mi ^ Trh- We find a 
small area of parameter space where a large enough baryon asymmetry is generated for Trh ~ 10^ GeV. 
It corresponds to large off-diagonal elements in to?, and therefore observable t ej. For Trh ^ 10^" 
GeV, leptogenesis can also work for smaller off-diagonal elements in to?. 
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Note added 

After this work was completed, related analyses appeared. 

-•^^If there are additional sources of flavour violation in the slepton masses, (e.g non-universal soft masses) this does not 
work. 
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7 Appendixrthe approximation and plots 

In this Appendix, the lightest eigenvalue and corresponding eigenvector of A4 are estimated, using an 
approximation borrowed from diagonalising neutrino mass matrices in R-parity violating theories. 

It is first convenient to scale some powers of the smallest Yukawa out of the hermitian matrix 

vtM-'^M'' = DyWL[m,]^VlDy''Vl[m,]VlDy' = A • (43) 
This can be written more compactly as 

4 = DY^AtDY-'ADY^ , (44) 

by defining 

A = Vl[m,]vl ^VlU*D,n^uWl = W*D„^^W^ , (45) 

where the matrix W — VlU is the rotation from the basis where the i//, masses are diagonal to the basis 
where the neutrino Yukawa matrix Yj^Yjy is diagonal. 
The matrix A can be written 

[Ak = (A.)-(A])-^(A.)fc(A*)fc , (46) 
fc 

where 

A,^^| y^Ayy, I . (47) 
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If the hierarchy in the yi is steeper than in the rrij, and/or that the angles in W are large, then 

|Aip»|A2|MA3p (48) 
SO the largest eigenvalue of A (= v^y\/\Mi\^) is 



2 2 2 2 2 2 



' '"vWF"|Aii| \Wi,m,\ ^^'^ 
with associated eingevector (normalised Ai): 

Ai ~ I yiAlJy2 I X (50) 
V yiA^i/ya I ^" 

In figure^] Mi is plotted as a function of 0^12 ~ logWi2 and ~ logWia, for j/i = lO^"', and using 
the central values of listed after eqn © . The precise definition is 

Wi2= COS Owi-i sin 9wi2 , Wi^ = sin Owii, 

with Owij = 10"i^7r/2 . (51) 

This Appendix contains many three dimensional plots of functions that will enter into the equation for 
the baryon asymmetry. The aim of these figures is to give a qualitative impression; quantitatively clearer 
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CO 13 



Figure 4: On the LHS, log[Afi/GeV] as a function of ~ loglWu] and ^13 ~ log[Wi3]. On the RHS, 
K = logio{rhi / m^) . Recall we need k ^ —1.4 to maximise the asymmetry. These plots are for central 
values of the neutrino masses: yi = 10^'', m2 = 8.2 x 10^'^ eV, and TO3 = 5.2 x 10^^ eV. 



contour plots of the baryon asymmetry are in the body of the paper. In figure01 there are three limiting 
values for Mi, corresponding to Mi ~ yfv'^/mi: Mi yfv'^/m^ when W13 1, Mi yfv'^/mi when 
Wi3, W12 0, and Mi ^ yfvl/m2 when W13 < TO2/TO3, W12 1. 

The fji decay rate can be evaluated with the eigenvector (|5()|l . which gives eqn H25(l . rhi has three 
limits — TOi,m2,m3 — depending ono the values of Wij. The logarithm of rhi/m^ is plotted on the RHS 
of figure^! rhi must be in the range given after eq. H16|l . which implies log (rhi/m^) ^ —1.4. 

Finally, the CP asymmetry e, eqn. (|20|) . can be evaluated with the eigenvector H5U|) to obtain 



3Ali ^ r [ADyA^DyA^]ii \ 3y2 ( EkW, 



where terms of order yi/y2 and yi/ys have been dropped. 5hmy oc ei/Mi is given in eqn (|26|l . and 
plotted on the LHS of figure[3 ei is plotted on the RHS; it peaks on the ridge of eqn H27(l because this 
is where the larger values of Mi and 5hmy overlap. 

The results in the remainder of the paper are based on the analytic approximations of this section. 
How reliable are these equations? The eqn (|44|l for A = yflAi/A"^]^^ is exact, but the formula we 
use for ei assumes hierarchical Mi, so it is consistent to assume this in solving for the eigenvalues and 
eigenvectors of A. The approximation is that the first column (or row) of yfOy^ ■ A ■ Dy^ is the lightest 
eigenvalue, multiplying its eigenvector. It breaks down if the elements of the second or third row/column 
become of order An, as one can see by writing the eigenvector in a basis rotated by a small angle from 
the eigenbasis. yiAi2/2/2, J/iAis/ys ~ An could occur if 
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Figure 5: On the LHS (RHS), Shmy (e) as a function of lji2 ~ ^05[T4^i2] and lui^ = loglWi^]. On the 
RHS, Ml is taken as a function of yi — 10""* and other inputs. Both plots are for central values of the 
neutrino masses. 

1. the Mi were of similar magnitude, rather than hierarchical. This is "unlikely", because the hier- 
archy in Dy is much steeper than in [rriu] ■ 

2. nil too small — if mi/m2, toi/tos < y|, then the terms being kept are smaller than the neglected 
ones. This is discussed in an appendix. 

8 Appendix: mi <C 1712/ 10 

In this paper, we assumed a hierarchical spectrum for the light neutrino masses: Ato^^^j = Wg, Am^^^ = 
m|, so the smallest neutrino mass mi is unlikely to be measured with anticipated data. However, it 
enters our formulae for the baryon asymmetry, as does the smallest Yukawa yi. In the body of the paper, 
we fixed mi = m2/10, and determined yi as a function of Mi and our weak scale parameters, by requiring 
Ml to be in the range allowed by leptogenesis. In this Appendix, we discuss the dependence of our results 
on mi. For most values of W12 and W13, mi is irrelevant because |Wiipmi <C |Wi2p'Ti2, jWispma. 
However, mi cannot be dropped from our analytic expressions, for W close to the identity. This is 
the area of parameter space where tjb is maximal; the remainder of the Appendix is restricted to this 
area of parameter space. We are interested in how 773 scales with mi, and in whether our analytic 
approximation is still valid. 

The maximum value of e, eqn H41|) . is independent of mi, if Mi is independent of mi. We have 
fixed Ml ~ Tj^jj, which determines yf as a function of Wi^m^ ~ mi. So varying mi allows j/i to vary: 
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mi ^ m2/100 would allow leptogenesis to work for yi ^ h^, which could be theoretically attractive. 

As discussed after eqn (|26|l . e approaches its upper bound (equivalently 5hmy ~ 1) on the peak in 
figure El where m^/mg ~ and < rnUrn^. As mi decreases, the area in W12, W13 space where 
Shmy 1 decreases, but the maximum value is unchanged. The numerical values of W12 and W13 where 
the maximum is reached will also decrease, making this parameter space increasingly "fine-tuned" {W 
very close to the identity is unlikely to be stable under renormalisation group running) . 

The i/n decay rate must have values in the range given after eqn (|16() . to ensure rj b b,s large as 
possible. We first concern ourselves with the upper bound: mi < 3 x 10^"^ eV. If Wf^m^ ~ Wi^ml and 
WijUij WiiTTii, as required to maximise Shmy, then from eqn H25|l . mi ~ Wi-^m^/mi. To maximise 
riB ^ Shmy /fni, requires W^^^ ~ m\/m\, so that the decay rate is slow enough, but Shmy is still 0(1). 
So as mi decreases from m2/10 to 10^'^m2, the area of the peak on the RHS of figure^will shrink, but 
the height is unchanged. 

For smaller values of mi, the asymmetry will decrease. This is because mi is small, so vb, production 
in the plasma is inefficient (see 

It is straightforward to check that our analytic approximation holds, in the shrinking area of pa- 
rameter space where W13 ~ mi/ms and W12 < mi/m2, provided that yf/mi <C yl/^^a- This is the 
condition that y^w^/mi is the lightest mass. So the analytic approximation fails as mi approaches 

Finally, the low energy prediction of the peak are independent of mi, because they follow from 
requiring that W 1. As mi decreases, W must approach the identity more and more closely, so 
becomes more precisely . But whether ~ t/^, or = U\ the expectation remains that t —> fij 
or r ^ 67 should be observable. 

So in summary, the magnitude of the baryon asymmetry on the peak of figure ^ is independent of 
mi for 10~^m2 < mi < m2/10. As mi decreases, the location of the peak shifts to smaller W13, and 
its area will shrink. We cannot say anything for mi < 10~'^m2: our analytic formulae indicate that rjB 
will decrease, but the approximation they are based on is unreliable. 
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